Biogenic crust dynamics on sand dunes 
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Sand dunes are often covered by vegetation and biogenic crusts. Despite their significant role 
in dune stabilization, biogenic crusts have rarely been considered in studies of dune dynamics. 
Using a simple model, we study the existence and stability ranges of difi'erent dune-cover states 
along gradients of rainfall and wind power. Two ranges of alternative stable states are identified: 
fixed crusted dunes and fixed vegetated dunes at low wind power, and fixed vegetated dunes and 
active dunes at high wind power. These results suggest a cross-over between two different forms of 
desertification. 



Sand dunes have been the subject of active research for 
many years, largely because of their fascinating shapes 
and dynamics \^]^. Current studies have increasingly 
addressed the question of sand-dune stability in relation 
to climate change and anthropogenic disturbances [B- 
Sand dunes are considered "stable" when they are 
fixed in place or are stationary 31|. Their stability is 
strongly affected by the degree of vegetation coverage. 
High coverage reduces the wind power at the dune sur- 
face and thereby acts to immobilize the dunes. The re- 
mobilization of fixed dunes, either by vegetation mortal- 
ity or clear-cutting, often has detrimental effects on the 
unique ecosystems that develop in stable dunes [1, 
leadiiig to alternative ecosystems associated with active 
sand [lOj . Active dunes may also pose a threat to human 
settlement as they can block roads and cover residential 
areas and agricultural fields [ll|, . 

Sand dunes are also stabilized by biogenic soil crusts. 
These crusts comprise a variety of organisms, including 
cyanobacteria, lichens and mosses, which live at the sur- 
face of desert soils Biogenic crusts enhance the ag- 
gregation of sand grains, prevent saltation, and reduce 
wind erosion. Since most sandy soils are located in dry- 
lands where the vegetation is patchy and generally sparse 
, the role of biogenic crusts in stabilizing dunes is im- 
portant and often crucial Despite their significance 
and vast presence in the Kalahari, Australian and Central 
Asia deserts lB"12|, soil crusts have rarely been consid- 
ered in studies of dune dynamics 18. [l9l|. 

Depending on wind power and precipitation level dif- 
ferent dune-cover states are observed. Figure [T] shows 
several typical states from regions of relatively weak 
winds. At very low precipitation levels (Fig. [1^), dunes 
are active due to low crust and vegetation coverage. As 
the precipitation increases, the dunes are gradually sta- 
bilized. At low precipitation levels, the dominant stabi- 
lizing agent is the soil crust (Fig. [T]d,c), while at high 
levels, the stabilizing agent is predominantly vegetation 
(Fig. [Hi). Although vegetation and biogenic crusts have 
similar effects in stabilizing dunes, they lead to ecosys- 
tems that differ vastly in their bioproductivity. 



Motivated by these observations, we ask whether the 
transition from crusted to vegetated dunes along the rain- 
fall gradient is gradual or abrupt, and how it is affected 
by the wind power. Studying these questions is signifi- 
cant for understanding desertification processes, i.e. pro- 
cesses involving the irreversible loss of vegetative biopro- 
ductivity (20j . To study these questions we introduce 
and analyze a new model, which extends an earlier model 
for vegetated dunes [111 to include crust dynamics. 




FIG. 1: Different types of sand dune cover along a rainfall 
gradient in Israel (p denoting the mean annual precipita- 
tion): (a) bare and active dune in hyper-arid region (Kasuy 
sands); (b) partially crusted and vegetated semi-active dune 
in arid region (Nizzana, 50% crust cover and 10% vegetation 
cover [IJl); (c) almost fully crusted and vegetated fixed dune 
in arid re gion (Halutza, 80% crust cover and 15% vegeta- 
tion cover ^21)' (f^) vegetated fixed dune in semi-arid region 
(Nitzanim); (e) the Israel-Egypt border region with crusted 
fixed dunes on the Israeli side alongside bare and mobile dunes 
on the Egyptian side. 
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The model consists of two coupled ordinary differential 
equations, 

2 

V = av{v + r]v)s - evDpVg{v)s - 'jD^v (la) 

~(j)yVb - ^yV 

b = ab{b + rib)s - ebDpbg{v)s (lb) 
~(f>bbv - ^ibb , 

for the fractions v and b of surface cover by vegetation 
and biogenic crust, respectively, where s — 1 — v — b>Q 
represents the remaining bare sand, and the over dot de- 
notes the time derivative. The first terms on the right 
sides of Eqs. (fTa|) and (flbl) represent logistic growth. Im- 
plicit in these growth forms is the assumption that the 
two life forms, crust and vegetation, locally exclude one 
another; the presence of crust in a given location prevents 
the germination of plant seeds, while the presence of veg- 
etation inhibits crust growth by blocking sunlight. The 
growth rates of vegetation, a^, and of biogenic crust, a;,, 
are assumed to have the following dependence on annual 
precipitation {p): 



ai{p) = 







(1 - e-^P-P'- 



P ^ Pi,min ; 



where i = v,b, Pmin is a precipitation threshold below 
which there is no growth, and ai^max is the asymptotic 
growth rate at high precipitation levels. This form is in 
accordance with field observations [2^. The parameters 
r]y and rjt represent spontaneous growth rates due to, for 
example, a bank of seeds and spores in the soil, respec- 
tively. 

Wind affects the vegetation and the crust both directly 
and indirectly. Indirect wind effects include sand trans- 
port, which leads to root exposure and burial of plants 
and crusts by sand. This process is represented by the 
second terms in ([1^) and (Ilt>)- The parameter Dp, the 
drift potential, is proportional to the potential bulk of 
sand that can be transported by the wind, and is given 
by Q: 



Dr, 



{U\U-Ut)), 



(3) 



where U is the wind speed and Ut is a threshold velocity 
that is necessary for sand transport (approximately 12 
knots for wind measured at 10 m above the ground). If 
U is measured in knots (1 knot « 0.5 m/s), the units 
of Dp are defined as vector units {VU). Dp can gen- 
erally be classified into low, intermediate and high en- 
ergy winds {Dp < 200VU, 200VU < Dp < 400VU, and 
Dp > 400VU, respectively 0). 

The function g{v) introduces a wind shielding effect 
created by vegetation. Observations indicate [2J| that 
when vegetative cover exceeds a certain value (vc), it in- 
duces a skimming flow in which sand is protected from 
direct wind action. This value depends on various prop- 
erties, such as plant shape and stem flexibility (25. [26|. 



Based on these studies, we chose a continuous step-like 
function for g{v): 
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g{v) = -(tanh(d(wc-t^)) + l). 



(4) 



such that g Q ioT V ^ Vc and g ^ \ ior v ^ Vc- The 
sharpness of g{v) is controlled by d. Since the indirect 
wind effect requires the availability of sand, the whole 
term is multiplied by s. 

Direct wind effects are restricted, in the model, to veg- 
etation and are represented by the third term in Eq. (jlap . 
This term accounts for stresses, such as increased evap- 
otranspiration and branch cutting. It does not have a 
parallel in Eq. (|lbp since crust can sustain very intense 
winds 27(1 . Wind drag is proportional to the square of 
the wind velocity, and therefore, this term is proportional 

i- n2/3 

to Dp . 

The parameters (j)y and 06 represent non-local in- 
teractions between vegetation and crust not accounted 
for by the first terms of the RHS. It is still debated 
whether these interactions are positive or negative. On 
one hand, crust supports vegetation growth as a result of 
the "source-sink" effect [l^: water runoff from the crust 
("source") flowing toward vegetation patches ("sink"). 
On the other hand, crust suppresses vegetation by pre- 
venting water infiltration in the vicinity of the roots of 
the plants. Biogenic crust is usually suppressed by plants 
due to litter from nearby plants which limits light and 
may destroy the crust if the litter is toxic [i^. Here, we 
assume that the negative relations are more significant. 
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Finally, /i^ and represent grazing stress 

The parameter values used in this study are based 
on Yizhaq et al. [3| for the equation of the vegeta- 
tion dynamics (Eq. (fTal) '). and on the studies of Bel- 
nap et al. 13|, l27| for the crust (Eq. 
merical values are: a^^max — 0.15 
50(mm-yr~^); c„ = 100 (mm ■ yr^^) ; ri^ — 0.2; 
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= (^yr ^j; ab^max 
cfc = 50 ■ 



= 0.015 (yr-i); 
Vb = 0.1; 



mm • yr " j; Cb — [mm ■ yr 
yr-i-VU-i); ^b - 0.01 (yr-i); ^lb = 
0.3; d= 15. 

The steady states of Eqs. (fTa|) and (|lbp and their sta- 
bility properties for low wind powers are presented in 
the bifurcation diagram shown in Fig. [2j The results 
are consistent with the general trend shown in Fig. [T] 
and reported in field observations [l3| : a low precipita- 
tion range (a) of bare active dunes; intermediate precipi- 
tation ranges (b,c) of dunes with mixed crust-vegetation 
coverage, semi-stabilized (b) or almost stabilized (c) ; and 
a high precipitation range (d) of stabilized vegetated 
dunes. In addition, the diagram predicts a bistability 
range (c) of vegetation-dominated dunes {v > b) and 
crust-dominated dunes (b > v). The bistability results 
from the negative vegetation-crust interactions assumed 
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in the model, which rehes on the conjecture that a crusted 
soil prevents the germination of plant seeds and also re- 
duces the infiltration of surface water into the soil, while 
a vegetated soil provides shading and possibly toxic ma- 
terials that inhibit the growth of crusts. 

A different type of bistability is known to exist in re- 
gions of strong winds and high precipitation [3l • This 
form of bistability results from the wind- shielding effect 
of the plants. The high wind power makes a bare dune 
active and prevents plant growth, despite the high pre- 
cipitation level. However, once the dune is vegetated, 
the resulting lower wind power allows its persistence. 
Thus, two forms of bistability, designated here as Type 
I and Type II, are possible. Type I is associated with 
the wind-shielding effect of vegetation and occurs at high 
precipitation and strong winds. In this case, the alterna- 
tive stable dune states are bare active dunes and vege- 
tated fixed dunes. Type II is associated with vegetation- 
crust competition and occurs at low precipitation and 
weak winds. Here, the alternative stable dune states are 
crust-dominated dunes and vegetation-dominated dunes. 
While bistability of Type I has been observed [l8| . 
observations of Type II have not yet been reported. 
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FIG. 2: Bifurcation diagram of dune-cover states at low wind 
powers. Solid and dashed lines represent, respectively, stable 
and unstable steady solutions of Eqs. (HJ with Dp — 120VU, 
obtained using a numerical continuation method. The red and 
blue lines represent, respectively, the vegetation-dominated 
and the crust-dominated states. The labels a,b,c,d refer to 
the corresponding panels in Fig. [T] They represent bare 
dunes (a); mixtures of crust and vegetation forming semi- 
stable dunes (b) and stable dunes (c); and stable vegetated 
dunes (d). 

FigureOshows the domains of the two bistability forms 
in the plane spanned by the precipitation p and the wind 
power Dp. The two domains are connected to form a 
continuous domain; proceeding from low to high p and 
Dp values, a cross-over from the bistability of Type II 
to Type I occurs. Bounding the continuous bistability 
domain are monostability domains of unvegetated dunes 



(bare or crusted) at low p or high Dp, and vegetated 
dunes at high p and low Dp. 
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FIG. 3: A phase diagram showing the monostability and 
bistability domains of dune-cover states. The shaded domain 
represents Type I bistability (active bare dune and fixed veg- 
etated dune) at high precipitation and wind power values, 
and Type II bistability (fixed vegetation-dominated dune and 
fixed crust-dominated dune) at low precipitation and wind 
power values. The difl^erent shades represent the difi^erence in 
sand-cover between the two stable states (calculated by sub- 
tracting the red line from the blue line in the bottom panel 
of Fig. [5)). The bistability domain is bounded by a monosta- 
bility domain of vegetated dunes at high precipitation, and a 
monostability domain of unvegetated dunes at low precipita- 
tion, crusted at low wind power and bare at high wind power. 

The existence of a biomass productive vegetation- 
dominated state and a less productive crust-dominated 
state, in the case of Type II bistability, implies the pos- 
sible occurrence of desertification, i.e. a state transition 
inducing bioproductivity loss, as well as the feasibility of 
rehabilitation of vegetation, a state transition resulting 
in bioproductivity gain. By "bioproductivity", we refer 
to the total amount of vegetative biomass. Such transi- 
tions can be triggered either by environmental variability, 
for example, precipitation or wind-power fluctuations, or 
by anthropogenic disturbances. The disturbance types 
that are necessary to trigger desertification or the reha- 
bilitation of vegetation can be determined by examining 
the positions of the two stable states in relation to the 
boundary between their basins of attraction, as Fig. 2] il- 
lustrates. Disturbances involving vegetation removal can 
induce desertification (transition from point B to point 
A) only at sufficiently low precipitation levels (panels (a) 
and (b)). At higher precipitation levels (panels (c) and 
(d)), the disturbance should also involve an increase in 
crust coverage (at the expense of sandy soil), a rather 
unlikely disturbance scenario. Rehabilitation of vegeta- 
tion (transition from point A to point B) at relatively 
low precipitation levels (panels (a) and (b)) cannot be 
triggered by crust removal only - planting is also neces- 
sary too. At higher precipitation levels (panels (c) and 
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(d)), crust removal alone can trigger such rehabilitation. 
These conclusions may be affected by other choices for 
the signs of the interaction coefficients 0„ and 0^. 

The desertification form discussed above should be dis- 
tinguished from that occurring in Type I bistability. In 
Type II bistability, both the productive and unproduc- 
tive states (i.e. vegetated and crusted) represent stable, 
immobilized dunes, while in Type I, the non-vegetated 
(unproductive) state represents a mobile dune. Thus, 
desertification in the case of Type I bistability not only 
involves the loss of vegetation (bioproductivity) but may 
also lead to detrimental effects associated with dune mo- 
bility. 
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FIG. 4: The phase space v, b for the bistability of Type 
II at increasing precipitation values [Dp — 120VU). The 
black dots A and B denote the crust-dominated state and 
the vegetation-dominated state, respectively. The line be- 
tween the dots represents the boundary between the basins 
of attraction of the two alternative stable states. The arrows 
represent disturbances or manipulations capable of inducing 
state transitions (see main text for details). 

The results of the model qualitatively agree with a phe- 
nomenon that has been observed at the Israeli-Egyptian 
border, where sand dunes on the Egyptian side are ac- 
tive, while on the Israeli side, dunes are semi-stabilized. 
This difference in dune activity is the result of a vast 
cover of biogenic crust on the Israeli side and its ab- 
sence on the Egyptian side. The phenomenon is clearly 
visible across the border line (see Fig. ^) due to dif- 
ferent albedo values for crust and sand. It was argued 
that biogenic crust is absent from the Egyptian dunes 
due to grazing activities that have led to crust tram- 
pling and erosion [30| . Analyzing this region in our 
model {p = 100 mm • yr~^), we find that grazing stress 
{Hv = ^ib = 0.01) yields bare dunes with low cover of 
crust and plants [h — 0.09, u = 0.11 — > s = 0.8), while 
absence of grazing {ji^ = /i;, = 0) yields crust-dominated 
dunes (6 = Q-b^v = 0.1 s = 0.4). Introducing grazing 
at low precipitation levels may, therefore, have a major 
effect on the fraction of bare dunes, doubling it in the 
numerical example presented above. 



Further analysis of the model shows that steady mixed 
dune-cover states may undergo oscillatory instabilities. 
We find that if non-local competition is assumed to be 
negligible ((/){, = 4>v =0), the oscillations emerge in a 
supercritical Hopf bifurcation at d = dc, where d is the 
parameter that controls the wind-shield steepness in Eq. 
(|4])). As d is increased beyond the Hopf bifurcation point, 
dc^ the amplitude of the oscillations increases and their 
frequency decreases. For d :s> dc (i.e. g{v) approaches a 
step function), the amplitude saturates and the dynamics 
take the form of relaxation-oscillations. 

We suggest the following mechanism for the oscilla- 
tions. Let us begin the description of a cycle with a state 
where v > Vc and b < Vc- Under such conditions, plants 
provide wind shield, which promotes the growth of both 
vegetation and crust at the expense of bare sand. As this 
growth proceeds, two processes become significant: the 
diminishing bare sand, which slows down the growth of 
both vegetation and crust; and the direct wind stress ef- 
fect that suppresses vegetation growth but does not affect 
the crust. The latter process favors the growth of crust 
at the expense of vegetation. As a result, the vegetation 
will reach a maximal value and start declining, while the 
crust will continue growing. The growth of crust and the 
decline of vegetation will continue until v < Vc- At this 
point, the sand drift effect will become dominant, and the 
crust will decline too, until there will be enough available 
sand for regrowth. Since the growth rate of vegetation 
is significantly larger than that of crust, vegetation will 
start growing at the expense of crust, until it exceeds Vc, 
and a new cycle will start over. 

In summary, a mathematical model was introduced to 
analyze the effect of biogenic crusts on sand dunes. Al- 
though simple, the model is able to capture important 
aspects of the complex dynamics of biogenic crusts and 
vegetation on sandy soils. Most significantly, it predicts 
a new form of bistability in which the two alternative 
stable states correspond to stabilized dunes with differ- 
ent proportions of vegetation and crust coverage. This 
bistability form (Type II) prevails at low precipitation 
and wind power values, and differs from the bistability 
of bare dunes and vegetated dunes at high precipitation 
and wind power values (Type I) The two bista- 

bility forms merge in the p — Dp parameter plane to form 
a single continuous domain with a small cross-over zone. 
The model sheds new light on the vulnerability of sandy 
regions to desertification and on the means to restore 
degraded vegetation. 
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